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ABSTRACT 



> 
<^ . 

, We present echelle spectrophotometry of the Galactic H ii region S 311. The data 

■rj" ' have been taken with the Very Large Telescope UVES echelle spectrograph in the 

\0 ■ 3100 to 10400 A range. We have measured the intensities of 263 emission lines, 178 

O ! are permitted lines of H^, 0° (deuterium), He", C°, C+, N^, N+, O", 0+, S+, Si", 

)Q- Si+, ArO and Fe"; some of them are produced by recombination and others mainly 

by fluorescence. Physical conditions have been derived using different continuum and 
r"| . line intensity ratios. We have derived He+, 0++ and 0+"'" ionic abundances from pure 

^ i' recombination lines as well as abundances from coUisionally excited lines for a large 

' I number of ions of different elements. We have obtained consistent estimations of ap- 

^ . plying different methods. We have found that the temperature fluctuations paradigm is 

consistent with the Te(He i) vs. rc(Hi) relation for H ii regions, in contrast with what 
^ [ has been found for planetary nebulae. We report the detection of deuterium Balmer 

lines up to Di5 in the blue wings of the hydrogen lines, which excitation mechanism 
seems to be continuum fluorescence. 

' Key words: ISM: abundances - H ii regions- ISM:individual: S 311 



1 INTRODUCTION 

S 311 — also known as N GC 2467 — is an H 1 1 region 
of the Sharpless catalogue iSharplessI Il959l) which is lo- 
cated outside the solar circle. This nebula forms part of 
the Puppis I association located at 4.0 kpc from the Sun 
jRusseil' 2003) and at a Galactocentric distance of 10.43 
(assuming a G alactocentric solar distance of 8.0 kpc). 
ert et all dlQSd) concluded that the radio morphology 
of S 311 is consistent with blister processes. There are few 
spectrophotometric studies of S 311 in the literature, most 
of them f orming part of analisys involving several H ll re- 
gions (^^gjHawlejlST^ JPeimbert. Torres-Peimbert fc Ravel 
Il978l: iKennicutt et aLlbOOOTl . All these works are based on 
the analysis of coUisionally excited lines (hereinafter CELs) . 

We have taken long-exposure high-spectral-resolution 
spectra with the Very Large Telescope (VLT) UVES echelle 



* Based on observations collected at the European Southern Ob- 
servatory, Chile, proposal number ESO 68.C-0149(A) 
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spectrograph to obtain accurate measurements of very faint 
permitted lines of heavy element ions in S 311. We have 
determined the physical conditions and the chemical abun- 
dances of S 311 with high accuracy. An important improve- 
ment in this work is the derivation of C"*"^ and O^"'' abun- 
dances from several pure recombination lines (hereinafter 
RLs) of C II and O ll, making use of high spectral resolution 
and avoiding the problems of line blending. 

Traditionally, the abundance studies for H ll regions 
have been based on determinations from CELs, whose emis- 
sivities are strongly dependent on the temperature varia- 
tions over the observed volume of the nebula. Alternatively, 
the emissivities of RLs are almost independent of such varia- 
tions and are, in principle, more precise indicators of the true 
chemical abundances of the nebula. The use of high resolu- 
tion echelle spectrographs has permitted our group to obtain 
deep high resolu tion spe ctra of bright Galactic H 1 1 regions 
fe.g.lEsteban et al. 199i ll999allbl F2004: Garcia- Roias et al, 
|200#7 and extragal actic H ll regions (e.g. .Esteban et alJ 
2002t |Peimbertll2003^ : all these works have found that abun- 
dance determinations from RLs are systematically larger 
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Table 1. Journal of observations. 



Telescope 


Date 


AA (A) 


Exp. time (s) 


8.2 m VLT 


2003/03/30 


3000-3900 


3x600 






3800-5000 


60,3x1800 


!J 




4750-6800 


3x600 


)' 




6700-10400 


60,3x1800 



than those obtained using CELs (the so called abundance 
discrepancy problem). One of the most probable causes of 
this abundance discrepancy is the presence of spatial varia- 
tions or fluctuations in the temperature structure of t he neb- 
ulae llTorres-Peimbert. Peimbert. fc Daltabuiti[l98ol) . Both 
phenomena can be related due to the different functional 
dependence of the line emissivities of CELs and RLs on the 
electron temperature, which is stronger — exponential — in 
the case of CELs. The temperature fluctuations have been 
parametrized traditionally by , the mean square tempera- 
ture fiuctuation of the gas ijPeimbertJIlQfiTl) . We have com- 
puted values from the comparison of abundances derived 
using CELs or RLs, and from the temperatures derived from 
CELs and recombination processes. 

The main aims of this work are to present the high- 
quality spectrophotometric data for S 311 obtained with 
the ESQ Very Large Telescope (VLT), to assess the old- 
fashioned abundance analysis of S 311 in the literature, cal- 
culate C^""" and O"*"^ abundances from recombination lines, 
and report the detection and measurement of weak deu- 
terium Balmer lines. 

In §§121 and 121 we describe the observations, the data re- 
duction procedure and the measurement and identification 
of the emission lines. In §2]we obtain temperatures and den- 
sities using several diagnostic ratios. In §|Klionic abundances 
are determined based on CELs, as well as on RLs. In §|Slwe 
discuss the results. Total abundances are determined in 
§ 13 In § (HI we discuss the detection of deuterium Balmer 
lines. Finally, in §§iniand 1101 we present the discussion and 
the conclusions, respectively. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were made on 2003 March 30 with 
the Ultraviolet Visu al Echelle Spectrograph, UVES 
JP'Odorico et al.ll200(il . at the VLT Kueyen Telescope in 
Cerro Paranal Observatory (Chile). We used the standard 
settings in both the red and blue arms of the spectrograph, 
covering the region from 3100 to 10400 A . The log of the 
observations is presented in Table Q 

The wavelength regions 5783-5830 A and 8540-8650 A 
were not observed due to a gap between the two CCDs used 
in the red arm. There are also five small gaps that were not 
observed, 9608-9612 A, 9761-9767 A, 9918-9927 A, 10080- 
10093 A and 10249-10264 A, because the five redmost or- 
ders did not fit completely within the CCD. We took long 
and short exposure spectra to check for possible saturation 
effects. 

The slit was oriented east- west and the atmospheric dis- 
persion corrector (ADC) was used to keep the same observed 
region within the slit regardless of the air mass value. The 
slit width was set to 3.0" and the slit length was set to 10" 
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Figure 1. Deep continuum-subtracted Ha image of S 311. The 
image was taken at the Nordic Optical Telescope (2.56m) in the 
Roque de los Muchachos Observatory, in La Palma (Spain) and is 
courtesy of A. R. Lopez-Sanchez. The black arrow indicates the 
position of the ionising star HD 64315. The white arrow indicates 
the position of the star NGC 2467-12 (see © The slit position 
is indicated as a white horizontal bar. 

in the blue arm and to 12" in the red arm; the slit width was 
chosen to maximize the S/N ratio of the emission lines and 
to maintain the required resolution to separate most of the 
weak lines needed for this project. The effective resolution 
at a given wavelength is approximately A A ~ A/8800. The 
centre of the slit was placed 126" south of the main ioniz- 
ing star HD 64315 (06e), covering the brightest region of 
S 311 (see Figure^. The reductions were made for an area 
of 3"x8.5". 

The spectra were reduced using the IRAF^ echelle re- 
duction package, following the standard procedure of bias 
subtraction, aperture extraction, flatfielding, wavelength 
calibration and flux calibration. The standard stars EG 247, 
C-32d9927 faamuv et alJll99l [1994) and HD 49798 were 
observed for flux calibration. 



3 LINE INTENSITIES AND REDDENING 
CORRECTION 

Line intensities were measured integrating all the flux in 
the line between two given limits and over a local contin- 
uum estimated by eye. In the cases of line blending, a mul- 
tiple Gaussian profile fit procedure was applied to obtain 
the line flux of each individual line. Most of these measure- 
ments were made with the SPLOT routine of the IRAF 
package. In some cases of very tight blends or blends with 
very bright telluric lines the analysis was performed via 

1 IRAF is distributed by NOAO, which is operated by AURA, 
under cooperative agreement with NSF. 
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Gaussian fittinR makinK u se of the Starlink DIPSO software 
l|Howard fc Murravl[l99ol) . 

Table |5| presents the emission line intensities of S 311. 
The first and fourth columns list the adopted laboratory 
wavelength, Ao, and the observed wavelength in the helio- 
centric framework, A. The second and third columns give 
the ion and the multiplet number, or series for each line. 
The fifth and sixth columns give the observed flux relative 
to H/3, F(\), and the flux corrected for reddening relative 
to H/3, /(A). The seventh colu mn gives the fractional e rror 
(Iff) in the line intensities (see lGarci'a-Roias et aljr2004l for 
details in error analysis). 

A total of 263 emission lines were measured; of them 
178 are permitted, 65 are forbidden and 19 are semiforbid- 
den (see Table I^J. Two Hi Paschen lines and one C 1 1 line 
are blended with telluric lines, making impossible their mea- 
surement. Several other lines were strongly affected by at- 
mospheric features in absorption, by internal reflections or 
charge transfer in the CCD, rendering their intensities un- 
reliable. Also, 23 lines are dubious identifications and one 
emission line could not be identified in any of the available 
references. Those lines are indicated in Table |5| 

The identification and adopted laboratory wavelengths 
of the lines were obtained follow ing several previous iden- 
tifications in the li terature (see ICarcfa-Roias et al.l 120041 : 
lEsteban et alj|2004l and references therein) 

We have assumed the standard extinct i on fo r the 
Milky Way (R„=3.1) parametrized bv ISeatonI (119791) . We 
have derived a logarithmic interstellar extinction coeffi- 
cient of c(_ff/3)=0.64 ± 0.04 dex was determined by fit- 
ting the observed /(H Balmer lines)/7(H/3) ratios (from 
H16 to H/3) and /(H Paschen lines) /J (H/3) (from P22 
to P7), to the theoretical ones computed by Storey & 
Hummer (1995) for Te = 10000 K and = 1000 cm~^ 
(see below). Hi lines affected by blends or atmospheric 
absorption were not considered. The derived value of 
c{HI3) is in very good agreeme nt with previous deter- 
minations in the same nebula: iHawlevI lll978f) derived 
c[H(5)—Q.Ql and 0.67 for two slit positions with offsets of 
96" north and 96" nort h, 35" west respectively from our slit 
positio n. Furthermore, IPeimbert. Torres-Peimbert fc Ravel 
i| 19781) derived c(J/"/3)=0.6 and 0.7 for slit positions 33" 
north, 12" west and 33" no rth, 106" west, r espectively. 
These last authors used the IWhitfordl il958l) extinction 
law, which is alm ost coincident with th e one we have 
adopted. Moreover. iKennicutt et all i2000t) . derived a value 
of c(H0) =Q.&7 using the average interstel l ar reddening 
curve from lCardelh. Clavton fc Mathij lll989h . lShaver et alJ 
il983 h observed two positions in S 311, one of them (position 
2), almost coincident with our slit position. They derived 
c{H (3)— 0.5 and 0.7 for their positions 1 and 2 respectively. 
We can conclude then that apparently there are no signifi- 
cant variations of the extinction inside S 311. 



4 PHYSICAL CONDITIONS 

The large number of emission lines identified and measured 
in the spectra allows us the derivation of physical condi- 
tions using different emission line ratios. The temperatures 
and densities are presented in Table |21 Most of the determi- 



Table 2. Observed and reddening-corrected line ratios 
[F(H/3)=100] and identifications. 



Aq •^obs err Notes 

(A) Ion Mult. (A) F(A)=' I{X)^ (%) 



3187.84 


He I 


3 


3188.45 


2.383 


4.486 


6 


3587.28 


He I 


32 


3588.08 


0.191 


0.294 


20 


3613.64 


He I 


6 


3614.45 


0.276 


0.422 


15 


3634.25 


He I 


28 


3635.08 


0.355 


0.537 


12 


3656.11 


H I 


H35 


3658.75 


0.070 


0.105 




3658.64 


H I 


H34 


3659.39 


0.103 


0.154 


34 


3659.42 


H I 


H33 


3660.26 


0.135 


0.202 


27 


3660.28 


H I 


H32 


3661.06 


0.180 


0.269 


21 


3661.22 


H I 


H31 


3662.10 


0.234 


0.350 


17 


3662.26 


H I 


H30 


3663.04 


0.177 


0.264 


22 


3663.40 


H I 


H29 


3664.25 


0.210 


0.314 


19 


3664.68 


H I 


H28 


3665.45 


0.255 


0.380 


16 


3666.10 


H I 


H27 


3666.91 


0.256 


0.382 


16 


3667.68 


H I 


H26 


3668.49 


0.270 


0.402 


15 


3669.47 


H I 


H25 


3670.26 


0.315 


0.470 


14 


3671.48 


H I 


H24 


3672.27 


0.375 


0.559 


12 


3673.76 


H I 


H23 


3674.57 


0.381 


0.567 


12 


3676.37 


H I 


H22 


3677.17 


0.496 


0.739 


10 


3679.36 


H I 


H21 


3680.16 


0.486 


0.722 


10 


3682.81 


H I 


H20 


3683.63 


0.545 


0.808 


9 


3686.83 


H I 


H19 


3687.63 


0.617 


0.914 


8 


3691.56 


H I 


H18 


3692.36 


0.777 


1.150 


7 


3697.15 


H I 


H17 


3697.97 


0.900 


1.329 


7 


3703.86 


H I 


H16 


3704.67 


1.024 


1.504 


6 


3705.04 


He I 


25 


3705.84 


0.407 


0.597 


11 


3711.97 


H I 


H15 


3712.79 


1.200 


1.756 


6 


3721.83 


[S HI] 


2F 


3722.68 


2.156 


3.145 


5 


3721.94 


H I 


H14 










3726.03 


[Oil] 


IF 


3726.88 


113.36 


165.15 


4 


3728.82 


[Oil] 


IF 


3729.63 


133.39 


194.14 


4 


3734.37 


H I 


H13 


3735.19 


1.812 


2.633 


5 


3750.15 


H I 


H12 


3750.97 


2.305 


3.330 


4 


3770.63 


H I 


Hll 


3771.47 


2.822 


4.048 


4 


3797.90 


H I 


HIO 


3798.74 


3.721 


5.289 


4 


3819.61 


He I 


20 


3820.47 


0.731 


1.031 


7 


3831.66 


S II 




3832.51 


0.026 


0.036 


29 


3833.57 


He I 


62 


3834.37 


0.117 


0.165 


10 


3835.39 


H I 


H9 


3836.23 


5.226 


7.338 


4 


3856.02 


Si II 


1 


3856.97 


0.046 


0.064 


20 


3856.13 


II 


12 










3867.48 


He I 


20 


3868.31 


0.084 


0.117 


12 


3868.75 


[Ne III] 


IF 


3869.62 


3.708 


5.152 


4 


3871.60 


He I 


60 


3872.67 


0.076 


0.105 


14 


3888.65 


He I 


2 


3889.50 


5.177 


7.149 


4 


3889.05 


H I 


H8 


3889.91 


7.991 


11.034 


4 


3920.68 


C II 


4 


3921.47 


0.034 


0.046 


25 


3926.53 


He I 


58 


3927.44 


0.098 


0.134 


11 


3964.73 


He I 


5 


3965.61 


0.668 


0.902 


4 


3967.46 


[Ne III] 


IF 


3968.34 


1.149 


1.549 


4 


3970.07 


H I 


H7 


3970.95 


11.680 


15.728 


4 


4008.36 


[Fe III] 


4F 


4009.25 


0.029 


0.039 


28 


4009.22 


He I 


55 


4010.14 


0.125 


0.167 


9 


4026.21 


He I 


18 


4027.10 


1.457 


1.927 


4 


4068.60 


[S II] 


IF 


4069.51 


1.710 


2.227 


4 


4076.35 


[S II] 


IF 


4077.25 


0.598 


0.777 


4 


4100.62 


D I 


D6 


4101.47 


0.025 


0.032 


30 


4101.74 


H I 


H6 


4102.64 


19.255 


24.792 


3 


4120.82 


He I 


16 


4121.75 


0.126 


0.161 


9 


4143.76 


He I 


53 


4144.68 


0.220 


0.280 


6 


4153.30 


II 


19 


4154.18 


0.019 


0.024 
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Table 2. continued 



Table 2. continued 



Ao 
(A) 




Ion 


Mult. 


Aobs 

(A) 




7(A)b 


err 

(%) 


4168, 


,97 


He I 


52 


4169.90 


0.035 


0.044 


24 


4169, 


,22 


O II 


19 










4267, 


,15 


C II 


6 


4268.11 


0.087 


0.108 


12 


4287, 


,40 


[Fe III 

L ^ J 


7F 


4288.35 


0.024 


0.029 


33 


4303, 


,61 


O II 


66 


4304.83 


0.030 


0.036 


27 


4303, 


,82 


O II 


53 










4339, 


,29 


D I 


D5 


4340.26 


0.036 


0.044 


22 


4340, 


,47 


H I 


H5 


4341.42 


38.644 


46.737 


3 


4359, 


,34 


fFe III 

L J 


7F 


4360.32 


0.022 


0.026 


35 


4363, 


,21 


\0 III] 


2F 


4364.17 


0.468 


0.562 


4 


4366, 


,89 


II 


2 


4367.95 


0.077 


0.092 


13 


4368, 


,15 


O I 


5 


4369.22 


0.032 


0.039 


26 


4368, 


,25 


O I 


5 










4387, 


,93 


He I 


51 


4388.90 


0.439 


0.522 


4 


4437, 


,55 


He I 


50 


4438.52 


0.063 


0.073 


15 


4471, 


,48 


He I 


14 


4472.49 


3.521 


4.056 


3 


4562, 


,60 


Mg I] 


1 


4563.63 


0.138 


0.153 


8 


4571, 


,10 


Mb I 


1 


4572.13 


0.120 


0.134 


9 


4630, 


,54 


N II 


5 


4631.65 


0.019 


0.021 


39 


4638, 


,86 


O II 


1 


4639.88 


0.026 


0.028 


30 


4641, 


,81 


O II 


1 


4642.85 


0.026 


0.028 


30 


4649, 


,13 


O II 


1 


4650.16 


0.023 


0.025 


33 


4650, 


,84 


O II 


1 


4651.81 


0.023 


0.025 


33 


4654, 


,12 


O I 


18 


4655.05 


0.013 


0.014 




4658, 


,10 


[Fc 111] 


3F 


4659.16 


0.184 


0.197 


7 


4661, 


,63 


O II 


1 


4662.60 


0.022 


0.024 


35 


4667, 


,01 


[Fc III] 


3F 


4667.91 


0.021 


0.022 




4701, 


,62 


[Fc III] 


3F 


4702.62 


0.050 


0.052 


18 


4711, 


,37 


Ar IV] 


IF 


4712.51 


0.012 


0.012 




4713, 


,14 


He I 


12 


4714.22 


0.394 


0.415 


4 


4754, 


,83 


[Fe III] 


3F 


4755.80 


0.041 


0.042 


21 


4769, 


,60 


[Fe III] 


3F 


4770.58 


0.019 


0.020 


36 


4777, 


,78 


[Fe III] 


3F 


4778.76 


0.015 


0.016 




4814, 


,55 


Fe III 


20F 


4815.64 


0.018 


0.018 




4815, 


,51 


S II 


9 


4816.67 


0.027 


0.027 


28 


4860, 


,03 


D I 


D4 


4861.08 


0.085 


0.085 


18 


4861, 


,33 


H I 


H4 


4862.40 


100.00 


100.00 


3 


4881, 


,00 


[Fc III] 


2F 


4882.13 


0.059 


0.059 


16 


4921, 


,93 


He I 


48 


4923.03 


1.117 


1.093 


4 


4924, 


,50 


[Fe III] 


2F 


4925.60 


0.020 


0.020 


38 


4924, 


,50 


O II 


28 










4930, 


,50 


[Fe III] 


IF 


4931.63 


0.017 


0.016 




4931, 


,32 


O III] 

J 


IF 


4932.28 


0.028 


0.027 


29 


4958, 


,91 


\0 III] 


IF 


4960.02 


44.315 


42.916 


3 


4985, 


,90 


[Fc III] 


2F 


4986.94 


0.142 


0.136 


12 


5006, 


,84 


o ml 


IF 


5007.96 


132.658 


126.110 


3 


5015, 


,68 


He I 


4 


5016.79 


2.328 


2.207 


3 


5035, 


,79 


[Fc II] 


4F 


5037.02 


0.023 


0.022 




5041, 


,03 


Si II 


5 


5042.25 


0.027 


0.025 




5047, 


,74 


He I 


47 


5048.92 


0.225 


0.211 


10 


5055, 


,98 


Si II 


5 


5057.12 


0.051 


0.048 


32 


5056, 


,31 


Si II 


5 










5191, 


,82 


[Ar III] 


IF 


5192.84 


0.056 


0.050 


30 


5197, 


,90 


[N I] 


IF 


5199.12 


0.389 


0.349 


7 


5200, 


,26 


[N I] 


IF 


5201.47 


0.330 


0.296 


8 


5270, 


,30 


[Fc III] 


IF 


5271.70 


0.124 


0.108 


16 


5517, 


,71 


[CI III] 


IF 


5518.91 


0.566 


0.461 


6 


5537, 


,88 


[CI III] 


IF 


5539.06 


0.437 


0.354 


7 


5577, 


,34 


[O I] 


3F 


5578.56 


0.072 


0.057 


25 


5754, 


,64 


[N II] 


3F 


5755.89 


1.315 


0.992 


4 



Ao 
(A) 



Ion 



Mult. 



Aobs 

(A) 



F{\Y /(A)*" 



err Notes 
(%) 



5875, 


,64 


He I 


11 


5876.96 


15.602 


11.343 


4 




5978, 


,83 


Si II 


4 


5980.31 


0.076 


0.054 


24 




6046, 


,44 


O I 


22 


6047.77 


0.034 


0.024 






6300, 


,30 


[OI] 


IF 


6301.75 


3.512 


2.307 


4 




6312, 


,10 


[S III] 


3F 


6313.48 


2.057 


1.348 


4 




6363, 


,78 


[O I] 


IF 


6365.24 


1.258 


0.815 


5 




6371, 


,36 


Si II 


2 


6372.85 


0.064 


0.041 


27 




6548, 


,03 


[N II] 


IF 


6549.54 


42.523 


26.484 


4 




6561, 


,04 


D I 


D3 


6562.43 


0.321 


0.200 


9 




6562, 


,82 


H I 


H3 


6564.26 


458.723 


284.845 


4 




6578, 


,05 


C II 


2 


6579.49 


0.196 


0.121 


12 




6583, 


,41 


[N II] 


IF 


6584.93 


131.646 


81.408 


4 




6678, 


,15 


He I 


46 


6679.64 


5.374 


3.263 


4 




6716, 


,47 


[S II] 


2F 


6717.97 


43.847 


26.434 


4 




6730, 


,85 


[S II] 


2F 


6732.35 


39.503 


23.752 


4 




6933, 


,91 


He I 


1/13 


6935.57 


0.017 


0.010 




g 


6989, 


,45 


He I 


1/12 


6991.12 


0.024 


0.014 






7002, 


,23 


I 


21 


7003.53 


0.179 


0.103 


9 


c 


7065, 


,28 


He I 


10 


7066.81 


3.780 


2.149 


5 




7092, 


,19 


CI] 




7093.57 


0.039 


0.022 


32 




7093, 


,97 


[Fe II] 




7095.54 


0.013 


0.007 




g 


7105, 


,85 


Si I 


70 


7107.13 


0.029 


0.016 




g 


7111, 


,47 


C I 




7113.06 


0.034 


0.019 


37 


g 


7135, 


,78 


[Ar III] 


IF 


7137.37 


17.966 


10.102 


5 




7155, 


,14 


[Fe II] 


14F 


7156.81 


0.029 


0.016 






7160, 


,58 


He I 


1/10 


7162.20 


0.048 


0.027 


27 




7231, 


,34 


C II 


3 


7232.89 


0.070 


0.039 


19 




7236, 


,19 


C II 


3 


— 


— 


— 


- 


c 


7281, 


,35 


He I 


45 


7282.98 


1.055 


0.580 


5 




7298, 


,05 


He I 


1/9 


7299.71 


0.039 


0.021 


33 




7318, 


,39 


[O II] 


2F 


7320.69 


1.985 


1.085 


4 




7319, 


,99 


[O II] 


2F 


7321.77 


6.024 


3.293 


5 




7329, 


,66 


[O II] 


2F 


7331.33 


3.311 


1.807 


5 




7330, 


,73 


[O 11] 


2F 


7332.41 


3.133 


1.710 


5 




7377, 


,83 


[Ni II] 


2F 


7379.55 


0.019 


0.010 






7423, 


,64 


N I 


3 


7425.40 


0.026 


0.014 






7442, 


,30 


N I 


3 


7444.02 


0.050 


0.027 


26 




7468, 


,31 


N I 


3 


7470.14 


0.156 


0.083 


10 




7499, 


,85 


He I 


1/8 


7501.53 


0.066 


0.035 


20 




7706, 


,74 


O I 


42 


7708.80 


0.036 


0.019 


35 




7751, 


,10 


[Ar III] 


2F 


7752.84 


4.854 


2.500 


5 




7782, 


,10 


Ca I] 




7784.06 


0.073 


0.038 


18 




7790, 


,60 


Ar I 




7792.83 


0.064 


0.033 


21 




7801, 


,56 


[Cr II] 




7803.60 


0.051 


0.026 


25 




7816, 


,13 


He I 


1/7 


7817.90 


0.106 


0.054 


14 




(oo I . 


■ OD 


[Co ij 




( OOU.DO 


U. iUO 


U.Uoo 


1 A 




7862, 


,75 


Fc II] 




7864.55 


0.023 


0.012 




g 


7875, 


,99 


[P II] 


iD-iS 


7877.68 


0.066 


0.033 


20 




7959, 


,70 


Co II] 




7961.62 


0.109 


0.055 


13 


g 


8046, 


,80 


Si I 


73 


8048.32 


0.048 


0.024 


27 


g 


8116 




He I 


4/16 


8116.93 


0.008 


0.004 






8150, 


,57 


Si I 


20 


8152.52 


0.077 


0.038 


19 




8184, 


,85 


N I 


2 


8186.78 


0.029 


0.014 






8188, 


,01 


N I 


2 


8189.93 


0.078 


0.038 


18 




8200, 


,36 


N I 


2 


8202.33 


0.017 


0.008 






8210, 


,72 


N I 


2 


8212.65 


0.036 


0.017 


35 




8216, 


,28 


N I 


2 


8218.25 


0.085 


0.042 


16 




8223, 


,14 


N I 


2 


8225.03 


0.095 


0.046 


15 




8243, 


,70 


H I 


P43 


8245.54 


0.068 


0.033 


20 




8245, 


,64 


H I 


P42 


8247.66 


0.094 


0.045 


15 
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Table 2. continued 



Table 2. continued 



Ao 
(A) 



err Notes 



Ion 



8247.73 
8249.97 
8252.40 
8255.02 
8257.85 
8260.93 
8264.28 
8266.40 
8267.94 
8271.93 
8276.31 
8281.12 
8286.43 
8292.31 
8298.83 
8306.11 
8314.26 
8323.42 
8333.78 
8334.75 
8345.55 
8359.00 
8374.48 
8387.77 
8392.40 
8395.98 
8413.32 
8437.96 
8446.25 
8446.36 
8446.76 
8451.00 
8459.32 
8467.25 



8502.48 
8518.04 
8528.99 
8665.02 
8680.28 
8683.40 
8686.15 
8703.25 
8711.70 
8718.83 
8727.13 
8733.43 
8750.47 
8776.77 
8845.38 
8850.62 
8862.79 
8888.71 
8893.87 
8914.77 
8946.05 
8996.99 
9014.91 
9019.14 



H I 
H I 
H I 
H I 
H I 
H I 
H I 
Ar I 
H I 
H I 
H I 
H I 
H I 
H I 
H I 
H I 
H I 
H I 
H I 
Fe II] 
H I 
H I 
H I 
Fe I 
H I 
Ca I] 
H I 
H I 
O I 
O I 
O I 
He I 
[Cr II] 
H I 
? 

He I 
H I 
He I 
He I 
H I 
N I 
N I 
N I 
N I 
N I 
N I 
[C I] 
He I 
H I 
He I 
He I 
Fe I] 
H I 
Fe I] 
VI] 
He I 
Fe II] 
He I 
H I 
Fe I] 







^ (AJ 








P41 


8249.60 


0.090 


0.044 


16 




P40 


8251.80 


0.081 


0.039 


17 




P39 


8254.35 


0.130 


0.063 


12 




P38 


8256.88 


0.104 


0.050 


14 




P37 


8259.72 


0.104 


0.050 


14 




P36 


8262.76 


0.112 


0.054 


13 




P35 


8266.18 


0.159 


0.077 


10 






8268.32 


0.092 


0.045 


15 




P34 


8269.78 


0.152 


0.074 


10 




P33 


8273.63 


0.086 


0.041 


16 


d 


P32 


— 


— 


— 




c 


P31 


8282.85 


0.270 


0.130 


8 


c 


P30 


— 


— 


— 


- 


c 


P29 


8294.22 


0.177 


0.085 


9 




P28 


8300.72 


0.252 


0.121 


8 




P27 


8307.94 


0.286 


0.138 


7 




P26 


8316.18 


0.301 


0.145 


7 




P25 


8325.28 


0.339 


0.162 


7 




P24 


8335.63 


0.434 


0.207 


6 






8336.68 


0.190 


0.091 


9 


g 


P23 


8347.42 


0.423 


0.202 


6 




P22 


8360.80 


0.629 


0.299 


6 




P21 


8376.33 


0.549 


0.260 


6 






8389.62 


0.095 


0.045 


15 


g 


P20 


8394.24 


0.634 


0.299 


6 






8397.84 


0.059 


0.028 


23 


c, 


P19 


8415.20 


0.713 


0.336 


6 


c 


P18 


8439.85 


0.818 


0.383 


6 




4 


8448.41 


0.823 


0.384 


6 




4 
4 


8448.93 


0.056 


0.026 


24 




6/17 


8453.21 


0.043 


0.020 


30 






8461.26 


0.193 


0.090 


9 




P17 


8469.15 


0.936 


0.435 


6 






8477.02 


0.032 


0.015 


39 




6/16 


8488.26 


0.034 


0.016 


37 




P16 


8504.33 


3.337 


1.538 


6 


c 


2/8 


8520.05 


0.076 


0.035 


18 




6/15 


8530.87 


0.060 


0.027 


22 


c 


P13 


8666.95 


2.136 


0.947 


6 




1 


8682.31 


0.064 


0.028 


21 




1 


8685.38 


0.152 


0.067 


11 




1 


8688.20 


0.044 


0.020 


29 




1 


8705.22 


0.051 


0.022 


26 


c 


1 


8713.75 


0.048 


0.021 


27 




1 


8720.90 


0.029 


0.013 








8729.19 


0.470 


0.205 


7 




6/12 


8735.41 


0.074 


0.032 


18 




P12 


8752.44 


2.672 


1.161 


6 




4/9 


8778.70 


0.074 


0.032 


18 




6/11 


8847.36 


0.100 


0.043 


15 






8852.74 


0.142 


0.060 


11 


g 


Pll 


8864.76 


3.597 


1.525 


6 






8890.93 


0.082 


0.035 


17 


g 




8895.90 


0.203 


0.086 


9 


c, 


2/7 


8916.58 


0.062 


0.026 


22 






8948.25 


0.088 


0.037 


16 


g 


6/10 


8998.95 


0.146 


0.060 


12 




PIO 


9016.97 


4.278 


1.764 


6 


d 




9021.38 


0.235 


0.097 


9 


g 



Ao 






A b 






err 


Notes 




Ion 


Mult. 


(A) 










9019.14 


Ca I] 














9029.07 


Fe I] 




9031.34 


0.145 


0.060 


11 


g 


9063.29 


He I 


4/8 


9065.18 


0.222 


0.091 


9 




9068.90 


[S HI] 


IF 


9070.92 


55.891 


22.855 


6 




9094.83 


C I 




9096.98 


0.139 


0.056 


12 


g 


9111.83 


Ca I] 




9113.95 


0.048 


0.020 


27 


g 


9113.60 


Fe I] 




9115.67 


0.108 


0.044 


14 


g 


9123.60 


[CI II] 




9125.72 


0.225 


0.091 


9 




9162.65 


Ni II] 




9164.72 


0.047 


0.019 


28 


g 


9210.28 


He I 


6/9 


9212.41 


0.189 


0.076 


10 




9226.62 


[Fe II] 




9228.58 


0.088 


0.036 


16 


g 


9229.01 


H I 


P9 


9231.11 


6.898 


2.770 


6 




9463.57 


He I 


1/5 


9465.71 


0.286 


0.113 


8 




9504.54 


CI] 




9506.97 


0.127 


0.050 


12 


g 


9o2d.10 


rle i 


6/8 


9525.46 


0.362 


0.143 


8 




9530.60 


[S HI] 


IF 


9533.12 


156.202 


61.749 


5 




9545.97 


H I 


P8 


9548.13 


7.956 


3.144 


7 


d 


9702.62 


He I 


75 


9705.23 


0.130 


0.051 


12 


c 


9824.13 


[CI] 


IF 


9826.49 


0.986 


0.389 


7 




9850.24 


[CI] 


IF 


9852.56 


2.766 


1.091 


7 




9876.67 


Fe I] 




9879.00 


0.233 


0.092 


9 


g 


10027.70 


He I 


6/7 


10029.98 


0.595 


0.235 


7 




10049.37 


H I 


P7 


10051.64 


17.935 


7.086 


7 




10286.70 


[S II] 


3F 


10288.95 


1.978 


0.782 


7 




10320.49 


[S II] 


3F 


10322.75 


2.235 


0.884 


7 




10336.41 


[S II] 


3F 


10338.96 


3.066 


1.212 


7 





" Where F is the observed flux in units of 100.00 = 
3.703 X 10-^3 ergs cm-^ s"!. 

^ Where / is the dereddened flux, in units of 100.00 = 
1.616 X 10~^^ ergs cm~^ s~^ and assuming C(H/3)=0.64 
dex. 

Affected by telluric emission lines. 
Affected by atmospheric absorption bands. 
^ Affected by internal reflections or charge transfer in the 
CCD. 

' Blend with an unknown line. 
8 Dubious identification. 



nations were carried out with the IRAF ta sk TEMDEN of 
the package NEBULAR JShaw fc DufouJ [l995'). 

The methodology followed for the derivation of rie 
and Te has been desc ribed in a previous paper (i.e. 
iGarcfa-Roias et al.ll200j) . In the case of electron densities, 
ratios of CELs of several ions have been used. The lastest 
version of NEBULAR (February 2004) uses the transi- 
tion probabilities recom mended by |Wiese_et_alJ ( 1996) and 
the collision strengths of iMcLauehlin .§^661111^)93.) for the 
[O ii] A3729/A3726 doublet ratio. These atomic data yield 
electron densities systematically lower than those deduced 
from the [S ll] A6716/A6731 dou blet ratio (see^teban et all 
2004: Garci'a-Roias ct alJ l2004l)^ Following t he arg uments of 
Copetti fc Writzl. (,200211 and IWang et all ll2004li we have 



adopted the transition probabilities from 




Zeippe 



n— I- ^^^^ 

and the collision strengths from IPradhari l)l97(fl for the 

[O ll] A3729/A3726 doublet ratio, which give electron den- 
sities that are in good agreement with the other density 
indicators. We have derived the [Fe ill] density from the in- 
tensity of the 6 brightest lines lines, which have errors less 
than 30 % and seem not to be affected by line blending, to- 



6 J. Garcia-Rojas et al. 



Table 3. Plasma Diagnostic. 



Parameter 



Line 



Value 



A^o {cm-3) [N i] (A5198)/(A5200) 
[O II] (A3726)/{A3729) 
[S II] (A6716)/{A6731) 
[Fe III] 

[CI III] (A5518)/(A5538) 
A^e (adopted) 

Tc (K) [N II] (A6548+A6583)/(A5755) 

[S II] (A6716+A6731)/(A4069+A4076) 
[O II] (A3726+A3729)/{A7320+A7330) 
To (low) 

[O III] (A4959+A5007)/(A4363) 
[Ar III] (A7136+A7751)/(A5192) 
[S III] (A9069+A9532)/(A6312) 
Tc (high) 
He I 

Balmer line/cont. 
Paschen line/cont. 



590+260 
260 ±110 

ogQ+140 

390 ± 220 

t:,;n+650 

310 ± 80 



9500 ± 250'' 
7200 +'^50 

(ZUU __gQ(j 

9800 ± eoo'' 
9550 ± 250 
9000 ± 200 
8800 t™« 
9300 ± 350 
9050 ± 200 
8750 ± 500 
9500 ± 900 
8700 ± 1100 



" Recombination contribution on the auroral lines has been taken 
into account (see text) 




8240 8260 

WavclcngLh (A) 



Figure 2. Section of the echelle spectrum including the Balmer 
(top) and the Paschen (bottom) limits (observed fluxes). 



Table 4. Atomic data used for some selected ions. 



Ion 



Coll. Strengths 



Trans. Probability 



0+ Pradhan (19761 

S+ Kocnan et al. (1996) 
S++ Taval fc Gupta (19991 



Zeippen (19821 
Keeimi^^al. (1993) 
Mendoza fc Zeippen (19821 



gether with the computations of iRodrfgued ll2002h . All the 
computed values of ric are consistent within the errors (see 
Table EJ. 

A weighted mean of n^iO ll), nc(Fe ill), no (CI III) and 
ne(S ll) has been used to derive T'o(N ll), re(0 ill), T'o(Ar III) 
and rc(S III), and iterated until convergence. So, for all the 
species we have adopted no=310 ± 80 cm~"^. We have ex- 
cluded no(N l) from the average because this ion is repre- 
sentative of the very outer part of the nebula, and probably 
does not coexist with most of the other ions. 

Electron temperatures have been derived from the ra- 
tio of CELs of several ions and making use of NEBULAR 
routines. In the cas e of ro ( S III) in stead of using the col- 
lision strengths of iCalavis. Mendoza fc ZeippenI lll995h . in- 
cl uded by default in NE BULAR, we have used the ones 
bv lTaval fc Guotal lll999h . The last set of collision strengths 
gives ro(S III) values more consistent with the rest of tem- 
perature determinations. Table0]shows the atomic data that 
we have changed in the last version of NEBULAR. 

To obtain To(0 ll) it is necessary to subtrac t the con- 
tribut ion to AA7320-I-7330 due to recombination; iLiu et all 
("2000^ find that the contribution to the intensities of the 
[O ll] AA 7319, 7320, 7331, and 7332 lines due to recombi- 
nation can be fitted in the range 0.5^T/10*^1.0 by: 



Jfl(7320 -f 7330) 
/(H/3) 



: 9.36 X {Tif-^'^ X 



H^ 



(1) 



where r4=r/10^. With this equation we estimate a contri- 
buti on of approxima tely 2% to the observed line intensities. 
iLiu et al] i2000ll also determined that the contribution 



to the intensity of the A 5755 [N ll] line due to recombination 
can be estimated from: 



/a (5755) 
/(H/3) 



= 3.19 X (T4) 



N^ 



H+ 



(2) 



in the range 0.5^ T/10'*^2.0. We have obtained a contri- 
bution of recombination of about 0.5%, that does not affect 
significantly the temperature determination. 

Figure |5| shows the spectral regions near the Balmer 
and the Paschen limits. The discontinuities can be easily 
appreciated. They are defined as Ic{Bac) — /c(A3646~) — 
7c(A3646+) and /o(Pac) = 7c(A8203-) - 7o(A8203+) re- 
spectively. The high spectral resolution of the spectra per- 
mits to measure the continuum emission in zones very 
near de discontinuity, minimizing the possible contamina- 
tion of other continuum contributions. We have obtained 
power-law fits to the relation between Ic{Bac) / I{Hn) or 
Ic{Pac) / I{Pn) and To for different n corresponding to dif- 
ferent observed lines of both series. The emissivities as a 
function of electron temperature for the nebular continuum 
and the H l Balmer and P asch en lines have been taken f rom 
iBrown fc Mathewj lll970ll and lStorev fc Hummeil re- 
spectively. The Ts{Bac) adopted is the average of the values 
using the lines from Ha to H 10 (the brightest ones). In 
the case of Te{Pac), the adopted value is the average of the 
individual temperatures obtained using the lines from P 7 
to P 13 (the brightest lines of the series), excluding P 8 
and P 10 because their intensity seems to be affected by sky 
abs orption. 

IPeimbert. Peimbert fc Luridianal i2002ll developed a 
method to derive the helium temperature. To (He i), 
in the presence of temperature fiuctuations. Assum- 
ing a 2-zone ionization scheme and the formulation of 
IPeimbert. Peimbert fc Luridianal (|20o3) we have derived 
To (He l)=8750 ± 500 K, which is highly consistent with 
To (Hi) assumed above. 

We have assumed a 2-zone ionization scheme for the 



Table 5. He+ abundance. 



Line He+ /H+ 



3819, 


.61 


775 


± 


54 


3888, 


.65 


759 


± 


30 


3964, 


.73 


825 


± 


33 


4026, 


.21 


818 


± 


33 


4387, 


.93 


837 


± 


34 


4471, 


.09 


792 


± 


24 


4713, 


.14 


821 


± 


33 


4921, 


.93 


800 


± 


32 


5875, 


.64 


780 


± 


31 


6678, 


.15 


793 


± 


32 


7065, 


.28 


770 


± 


38 


7281, 


.35 


835 


± 


42 



Adopted 795 ± 9*= 



In units of lO"", for T3889=2.52 ± 0.44, and (2=0.034 
± 0.010. Uncertainties correspond to line intensity errors. 

^ It includes all the relevant uncertainties in emission 
line intensities, no, t^ss9 ^'Hd ■ 

derivation of ionic abundances (see §01. We have adopted 
the average of electron temperatures obtained from [N ll] 
and [O ll] lines as representative for the low ionization zone, 
and the average of the values obtained from [O ill], [Ar ill] 
and [S III] lines for the high ionization zone (see Table I^J- 

5 IONIC ABUNDANCES 

5.1 He^ abundance 

We have measured 50 He I emission lines identified in 
our spectra. These lines arise mainly from recombination 
but they can be affected by coUisional excitation and 
self-absorption effects. We have determined the He^/H"'' 
ratio using the e ffectiv e recombination c oefficients o f 
IStorev fc Hummed ^1995^ for Hi and th ose of lSmitj (Il996l) 
and iBei^Mnin ^^ilhnaii & Smith^ fm^ for H e I. The coUi- 
sional contribution was estimated from Sawov fc BerringtonI 
l| 19931) and iKingdon fc Ferlandl lll995l) . and the optical 
depths i n the triplet lines were derived from the computa- 
tions bv lBeniamin. Skillman fc Smit^ (|200^. From a max- 
imum likelihood method fe.g. IPeimbe^Peiinbert fc Ruid 
\200(t . using ne=310 ± 80 cm"^ and T(0 ll-hiii)=9600 ± 
450 K (see §|nil, we have obtained He+/H+=0.0795 ± 0.0009, 
T3S89=2.52 ± 0.44, and t^=0.034 ± 0.010. In Tablel^lwe have 
included the He^/H^ ratios we have obtained for the indi- 
vidual He I lines not affected by line blending and with the 
highest signal-to-noise ratio. We have excluded He I A 5015 
for the same reasons outlined bv lEsteban et afl i2004ft . We 
have done a optimization of the values in the table, and 
we have obtained a parameter of 8.15, which indicates a 
reasonable goodness of the fit for a system with nine degrees 
of freedom. 

5.2 Ionic Abundances from CELs 

Ionic abundances of N+, 0+, 0++, Ne++, S+, S++, C1+, 
Cl"*"^, Ar"'"^ and Ar^+ have been determined from CELs, 
using the IRAF package NEBULAR (except for Cl^, see 
iGarcia-Roias et alJl2004|) . Additionally, we have determined 
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Table 6. Ionic abundances from coUisionally excited 
lines''. 



Ion 


(2=0.000 


4^=0.038 ± 0.007 


N" 


c; 74 -1- n nfi 




N+ 


7.25 ± 0.05 


7.38 ± 0.06 


QO 


6.74 ± 0.06 


6.87 ± 0.06 


0+ 


8.26 ± 0.07 


8.40 ± 0.08 


0++ 


7.81 ± 0.04 


8.05 ± 0.06 


Ne++ 


6.81 ± 0.05 


7.07 ± 0.07 


S+ 


6.03 ± 0.05 


6.15 ± 0.06 


S++ 


6.68 ± 0.07 


6.95 ± 0.09 


C1+ 


4.56 ± 0.05 


4.67 ± 0.05 


C1++ 


4.85 ± 0.05 


5.08 ± 0.05 


Ar++ 


6.08 ± 0.04 


6.28 ± 0.06 


Ar3+ 


3 42+0-18 




Fe++ 


5.05 ± 0.06 


5.30 ± 0.08 



In units of 12+log(X'"/H+). 

the ionic abundan ces of Fe^"*" followi n g the methods and 
data discussed in iGarcia-Roias et al.l (|200^. Ionic abun- 
dances are listed in Table |H| and correspond to the mean 
value of the abundances derived from all the individual lines 
of each ion observed (weighted by their relative intensities) . 

To derive the abundances for — 0.038 (see § |SJl 
we us ed the a bundances for ^'^=0.00 and the formula- 
tion of lPeimber t (1967) and lPeimbert fc Costerol ^1969^ for 
>0.00. To derive abundances for other values it is pos- 
sible to interpolate or extrapolate the values presented in 
Table H 

Many [Fe ll] lines have been identified in our spectra, 
but all of them are severely affected by fluorescence effects 
jRodn'guedll999l: IVerner et al.ll200dl . The only [Fe ll] line in 
the spectral range 3100 A to 10400 A which is not affected 
by fluorescence effects is the [Fe ll] A8617 A line, but un- 
fortunately it is in one of our observational gaps. We have 
measured [Fe ll] A715 5, a line which is not much affected 
by ffuorescence effects JVerner et alJl2000H . but it has a high 
observational error. Therefore, it was no possible to derive 
a reliable value of the Fe''"/H"*' ratio. The calculations for 
Fe''""^ have been done wit h a 34 l e vel m odel-atom that uses 
the collision st rengths of IZhand (Il996l) and the transition 
probabilities of lOuine^ lll996ll . We have used 6 [Fe ill] lines 
that do not seem to be affected by line-blending and with 
errors less than 30 %. We find a mean value and a standard 
deviation of Fe++/H+ = (1.115 ± 0.153) x 10"''. Adding er- 
rors in To and Ue we finally obtain 12-|-log(Fe^"'"/H^)=5.05 
± 0.06. The value of the Fe++ abundance for t^ >0.00 is 
also shown in Table 15] 

5.3 Ionic Abundances from Recombination Lines 

We have measured a large number of permitted lines of 
heavy element ions such as O I, O ll, C I, C ll, S ll, N I, N ll, 
Ar I, Si I, Si ii, and Fe i most of them detected for the first 
time in S 311. Those permitted lines produced by recombi- 
nation can give accurate determinations of ionic abundances 
because their intensities depend weakly on electron temper- 
ature and density. Unfortunately most of the permitted lines 
are affected by fluorescence effects or are blended with tel- 
luric emiss i on lin es making their intensities unreliable; also 
iRuiz et al.l (|200^ have shown that, to determine the abun- 
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Table 7. C"'"+/H+ abundance ratio from C II lines 



Table 9. parameter 







/(A)//(H/3) 


C++/H+ 


(xlO-^) 


Mult. 


Ao 


(xlO-2) 


A 


B 


2 


6578.05^ 


0.121 ± 0.015 


129 ± 15 


23 ± 3 


3 


7231.12 


0.038 ± 0.007 


1037 ± 197 


15 ± 3 


4 


3920.68 


0.046 ± 0.012 


400 ± 100 


130 ± 33 


6 


4267.26 


0.108 ± 0.013 


10 ± 1 


10 ± 1 



Method 




t2 




0++ (R,/C) 


0.040 


± 


0.008 


He+ 


0.034 


± 


0.010 


Bac-FL 


0.002 


± 


0.022 


Pac-FL 


0.019 


± 


0.026 



Adopted 0.038 ± 0.007 



Adopted 



10 ± 1 



Affected by a telluric emission line 



Table 8. 0++/H+ ratio from O II permitted lines^ 







/(A)//{H/3) 


0+- 


7h+ (Xl0-5)b 


Mult. 


Ao 


(xlO-2) 


A 


B C 


1 


4638.85 


0.028 ± 0.008 


14 ± 4 


14 ± 4 




4641.81 


0.028 ± 0.008 


13 ± 4 


12 ± 4 




4649.14 


0.025 ± 0.008 


13 ± 4 


12 ± 4 




4650.84 


0.025 ± 0.008 


11 ± 4 


11 ± 4 




4661.64 


0.024 ± 0.008 


10 ± 4 


10 ± 4 




Sum 




12 ± 1 


12 ± 1 


2 


4366.89 


0.092 ± 0.012 


25 ± 3 


18 ± 2 


19 


4153.30 


0.024: 


794:/31: 


28:/28: 28:/29: 




Adopted 






12 ± 1 



Only lines with intensity uncertainties lower than 40 % have been 
considered 

^ Abundances of multiplet 1 have been corrected from NLTE effects 
(see text) 

dances, it is important to measure all the lines of a multiplet, 
because for low densities there could be an anomalous distri- 
bition of the line intensities within the multiplet. Detailed 
discussions on t he mechanism of formation of the permit- 
ted lines are in lEsteban et alJ il998l l2004l and references 
therein) . 

We have been able of measuring ionic abundance ra- 
tios of 0'''^/H^ and C^'''/H^ from pure recombination 

II and C ll lines respectively, from multiplet 1 for O ll 
(see Fip^ure ^ and from multiplet 6 for C ll [see Figure 

1 of lEsteban et all (|200^]. We have computed the abun- 




4640 4650 4660 4670 

Wavelength (A) 

Figure 3. Section of the echelle spectrum showing the lines of 
multiplet 1 of O II (observed fluxes). 



dances for ro(High)=9050 K and nc=310 cm ^. Atomic data 
and methodology for the C abundance are the same as in 
iGarcfa-Roias et al.l J2004) . For the Multiplet 1 of O II it has 
been shown that for densities Ue < 10000 cm~^ the upper 
levels of the transition are not in LTE, and if one uses only 
one line to determine th e abundances one can have errors as 
large as a factor of 4 llR.uiz et |2(^(^ j |) ; instead we used 
the p rescription presented bv IPeimbertT Peimbert fc Rui'3 
1I2OO3) to calculate those populations; these abundances 
show very good agreement among themselves and with the 
abundance determined using the sum of all the lines, "Sum", 
that is not expected to be affected by this effect. Tables Q 
and |H| show the abundance ratios. 



6 TEMPERATURE VARIATIONS 

It is well known that under the assumption of a con- 
stant temperature, RLs of heavy element ions yield 



in H II regions (e.g. IPeimbert. Storev & Torres-Peimberl 


ngga lEsteban et alj Il99d: lEstebanl l2002l 


lEsteban et al. 


l200llGarci'a-Roias et al."2004": 'Tsamis et al.' 


'200:?, and ref- 



(1980) proposed the presence of spatial temperature fluc- 
tuations (parametrized by t'^) as the cause of this dis- 
crepancy, because CELs and RLs emissivities have dif- 
ferent depende nces on the electron temperature. On the 
other hand, iPeimberd (^^) proposed that there is a di- 
chotomy between Te derived from the [O ill] lines and 
from the hydrogen recombination continuum discontinu- 
ities, which is strongly correlated with the discrepancy be- 
tween GEL and RL abundances, so the comparison be- 
tween electron temperatures obtained from both methods 
is an additional indicator oi . A complete form ulation of 
tempe rature fluctuations has been devel oped bv [ P eimberti 
Jl967f) . [Peimbcrt & Gostcro (Il969l) and iPeimberti Jl97ll) . 
We have assumed a two-zone io nization scheme, and we have 
follow ed th e re-formulation of [Pe imbert^ Peirnbe rt fc Rui3 
J2000l) and iPeimbert. Peimbert fc Luridianal ^2002^ to de- 
rive the value of comparing Tc{Bac) and Tc{Pac) with 
the combination of rc([0 11]) and re([0 ill]), t'^(Bac - FL) 
and t^jPac — FL), respectively , usin g equation (Al) of 
IPeimbert. Peimbert fc Luridianal (l2002h . Table El shows the 
different values obtained as well as the adopted value, 
t'^=0.038 ± 0.007 which is the weighted average of O^""" and 
He^ values which are rather consistent and show the lowest 
uncertainties. 

On the other hand, as it can be seen from Table (51 
t'^{Bac - FL) and t'^{Pac - FL) are signiflcantly lower 
than the other values. One possible explanation is that 



Echelle Spectrophotometry of S 311 9 



Table 10. Continuum determinations'' 



A(A) log (j(A))//(H/3) 





Atomic 


Observed 


Scattered light 


3640 


-2.235 


-2.106 ± 0.009 


-2.699 ± 0.033 


3670 


-3.057 


-2.469 ± 0.025 


-2.598 ± 0.033 


4110 


-3.185 


-2.577 ± 0.024 


-2.700 ± 0.031 


4350 


-3.213 


-2.676 ± 0.022 


-2.825 ± 0.031 


4850 


-3.245 


-2.767 ± 0.015 


-2.943 ± 0.022 


6650 


-3.310 


-2.870 ± 0.013 


-3.066 ± 0.020 


8190 


-3.354 


-3.059 ± 0.004 


-3.367 ± 0.009 


8260 


-3.883 


-3.262 ± 0.003 


-3.381 ± 0.004 



in units of (A~^) 



the nebular continuum could be affected by dust scat- 
tered light. To explore that possibility we have derived the 
atomic continua contribution, which includes the free-free 
and free-bound continua of the hydrogen and helium atoms 
a nd the two-quantum conti nuum, from the computations 
of iBrown fc Mathewj (|l97nl ^ for T^=762Q K, no=310 cm"^ 
and He^/H"'"=0.0795 (see section IS!T]l . The temperature of 
7620 K we have assumed is that which implies t^{Bac~ FL) 
and t^{Pac — FL) equal to the final adopted value. Ta- 
ble 1101 shows the observed and the expected atomic con- 
tinua, and the derived scattered light contribution. From 
these data it is easy to derive the contribution of dust scat- 
tered light to the continuum near the Balmer lines and 
the Balmer and Paschen jumps. Including observational un- 
certainties, we have derived that a contribution between 
10% and 30% to the Balmer jump of the scattered contin- 
uum integrated light is enough to explain the high Tc{Bac) 
obtained from the discontinuity and, therefore, the lower 
t^{Bac — FL). The shape of the spectral distribution of 
the scattered light is then similar to a B3 V star. Assum- 
ing that the ionising star of S 311 -HD 64315- is a main 
sequence 06e, we have modelled the optical flux distribu- 
tion using FASTWIND, a n spherically symmetric, NLTE 
model atmosphere cod e l|Santolava-R,ev. Puis, fc Herrerol 
Il997l: IPuls et alJ l200!^^ and the correspondi n g Teg and 
logg derived from iMartins. Schaerer fc HillieJ ll2005h cali- 
brations. From this model, we have obtained that the stel- 
lar Balmer jump is about 10% of the stellar continuum. 
Anyway, our slit position is far from HD 64315 (see Fig- 
ure and other late O and early B stars may contribute 
to the continuum scattered light, so the contribution of 
the discontinuities would be even higher, and therefore, the 
nebular temperature determ ination would decrease. In fact, 
iFeinstein. fc Vazaued lll98Sl) list a B8.5 star near our slit 
position (their star labelled as NGC 2467-12), whose posi- 
tion is also indicated in Figure 1. It can be seen that, as 
expected, the continuum scattered light increases monoton- 
ically towards the blue. A more detailed study of the proper- 
ties of the dust in S311 (albedo, reddening and geometrical 
distribution) is needed to solve this problem in an appropi- 
ate way, but t his is outside the sc ope of this paper. 

R.ecentlv. IZhang et alJ ('2005) have computed To (He l) 
for 48 planetary nebulae from He I recombination line ratios 
(He IA7281/A6678). These authors have found that tempera- 
ture fluctuations do not predict the general behaviour of the 
TefHe l) vs. re(Hl) diagram. Using the expression given by 
IZhang et aO J2005li we have derived rc(He i)=7960 ± 1000 




5000 10* 1,5x10' 2x10' 

r/IIo I X7281/A6678) (K) 



F igure 4. TeCRe l) versus Te(Hl). This figure is like Figure 5 
of lZhang et al.l l2005^ for PNe (open squares), including data for 
H o regions (filled circles). The solid diagonal line is a y = x plot, 
and the dotted lines show the variations of Te(Hl) as a function 
of Te(He l) for diff'erent values. The cross on the lower right 
part of the diagram represent typical uncertainties for PNe data. 



K, for S 311, which is lower, but consistent within the errors 
with our Te (He l) and Te(Hl). I n Figure|31we have compared 
results from lZhang et al"! (l2005l) for planetary nebulae (PNe) 
with our results for H ll regio ns, in which we have incl uded 
S311 (this paper), NGC 3576 JOarcfa-Roias et al.l20o3) and 
still unpublished echelle VLT spectra of M8, M17, M20, M16 
and NGC 3603. It can be seen that most of the H ii regions, 
exce pt for the Orion nebula (lEsteban et alj|2004l) and 30 
Dor iPeimbertil2003) . do not follow the behaviour of the 
bulk of PNe and do not contradict t he temperature fluctu- 
ations paradigm. IZhang et al.l (j200^ solve the problem of 
the anomalous position of PNe in the Te(He l) vs. re(Hl) 
diagram proposing the presence of a small amount of H- de- 
flcient material for most of the sample nebulae, in the con- 
text of the scenario of H- deflcient inclusions constructed by 
iPeauignot et alJ i2003l) . That scenario seems to explain suc- 
cessfully the large abundance discrepancies reported in some 
PNe. However, Figure |1] suggests that the behaviour of H ll 
regions is qualitatively different to that of PNe and that the 
position of most H ll regions in the diagram is consistent 
with the classical paradigm. T his result is in agreement 
with the suggestions outlined bv lEstebanl i2002h who pro- 
pose that the processes producing abundance discrepancies 
in H II regions and PNe -at least the extreme cases of PNe- 
could be different. In the case of H ii regions the observa- 
tional evidences are still consistent with the scheme, while 
in the case of the extreme PNe this scheme seems to fail. 
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Table 11. Total Gaseous Abundances. 



Element 


S 311 (this work) 
(2=0.000 (2=0.038 ± 0.007 


Orion^ 
(2=0.022 ± 0.002 


Sunl^ 


S 311-Orion 


He 


10.99 ± 0.02 


10.97 ± 0.02 


10.988 ± 0.003 


10.98 ± 0.02 


-0.018 


C 


8.38 ± 0.05 


8.38 ± 0.07 


8.42 ± 0.02 


8.39 ± 0.05 


-0.04 


N 


7.43 ± 0.06 


7.61 ± 0.07 


7.73 ± 0.09 


7.78 ± 0.06 


-0.12 


O 


8.39 ± 0.05 


8.56 ± 0.06 


8.67 ± 0.04 


8.66 ± 0.05 


-0.11 


o° 


8.54 ± 0.10 


8.57 ± 0.05 


8.65 ± 0.03 




-0.08 


Ne 


7.79 ± 0.13 


7.98 ± 0.14 


8.05 ± 0.07 


7.84 ± 0.06 


-0.07 


S 


6.77 ± 0.06 


7.02 ± 0.08 


7.22 ± 0.04 


7.14 ± 0.05 


-0.20 


CI 


5.03 ± 0.06 


5.22 ± 0.07 


5.28 ± 0.04 


5.23 ± 0.06 


-0.06 


Ar 


6.43 ± 0.07 


6.56 ± 0.08 


6.62 ± 0.05 


6.18 ± 0.08 


-0.06 


Fe 


5.17 ± 0.11 


5.44 ± 0.13 


6.23 ± 0.08 


7.45 ± 0.05 


-0.79 



" Gas abundances from lEsteban et alJ J2223) 
IChristensen-DalsgaardI il998^ : IXsplund. Grevesse fc Sauvall <200,'J) 
0++/H+ from RLs and 0+/H+ from CELs and t^ 



7 TOTAL ABUNDANCES 

We have adopted a set of ionization correction factors 
(ICFs) to correct for the unseen ionization stages and then 
derive the total gaseous abundances of the elements we 
have studie d. We have a dopted the ICF scheme used by 
[Gar cfa-Roia s et alJ i2004l) for all the elements except for 
carbon, neon, chlorine and iron. 

The absence of He ll lines in our spectra indi- 
cates that He"''^/H^ is negligible. However, the total he- 
lium abundance has to be corrected for the presence 
of neutral helium. Based on th e ICF (He'') given by 
IPeimbert. Torres-Peimbert fc Rui3 lll992ll . and our data 
ICF(He°) amounts to 1.22 ± 0.05 for = 0.00 and 1.16 
± 0.04 for > 0.00. 

We have derived the 0/H ratio both from CELs and 
from the combination of 0"'"^/H^ ratio from RLs and 
0~'"/H''' ratio from CELs and the assumed t^. 

For carbon we hav e adopted the ICF de rived from pho- 
toionization models bv lCarnett et all il9&. " 



For neon the ICF proposed bv lPeimbert fc Costerol 
1969i ) given by: 



N{Ne) 

Apr 



N{0+) + N{0++)\ N{Ne++) 



iV(0^ 



iV(H+) 



(3) 



has been generally used. Nevertheless this ICF un- 
derestimates the Ne/H abundance for nebulae of low 
degree of ionization becaus e a con siderabl e fraction of 
Ne"*" co exists with O"*"^ (see iTorres-Pcimbcr t. fc Peimberd 
Il977l: IPeimbert. Torres-Peimbert fc Ruij Il992h . For 
S 311 based on the O^/ O rati o and the data by 
iTorres-Peimbert. fc Peimberd (Il977l) . we estimate that 
the ICF(Ne) should be 0.4 ± 0.1 dex higher than that 
provided by the previous equation. 

We have measured lines of two ionization stages of chlo- 
rine: Cl^ and Cl"*"^. The CI abundance has been assumed to 
be equal to the sum of these ionic abundances without taking 
into account Cl^"*" fraction. This assumption seems reason- 
able taking into account the small Cl^^ / Cl^"*" ratio found for 
M17 (~0.0 3. seelEsteban et'aI]ll999a^ . for the Orion nebula 
f^ 0.04. seelEsteban et al.ll2004h . and for NGC3576 (~0.02, 
see lGarci'a-^oia^^n]20M) . and the lower ionization de- 
gree of S 311 with respect to those nebulae. 



We have measured lines of two stages of ionization of 
iron: Fe^ and Fe"""^, but in ? l5.2lwe have shown that F e ^/H"* " 
ratio is not reliable. Recently, iRodrfgu'ez' fc RubinI diooi) 
have derived an ICF from a least-squares fit to the results of 
a set of models in which it is represented x(0''")/x(Fe^^)) 
the ratio of ionization fractions of and Fe^"'', respec- 
tively, as a function of the degree of ionization given by 
0"''/0'''^. The ionization correction scheme they have de- 
rived is as follows: 



iV(Fe) 



0.9 



N{0+) 
N{0++) 



jV(Fe)++ N{0) 
N{0)+ ^ N{R)' 



(4) 



In Table im we show the total abundances obtained for 
our slit position in S 311 for (^=0.00 and (2=0.038 ± 0.007. 



8 DEUTERIUM BALMER LINES 

We have detected, for the first time in an H ll region outside 
the solar circle, the four brightest deuterium Balmer lines. 
Da, D/3, D7 and DS as very weak lines in the blue wings 
of the corresponding H I Balmer lines (see Figure |3 . The 
apparent shift in radial velocity of these weak lines with 
respect to the hydrogen ones is -82.7 km s"^ (see TablelT^. 
which is in excellent agreement with the isotopic shift of 
deuterium, —81.6 km s~^. We have discarded these weak 
features as high velocity components of hydrogen for the 
following reasons: 

• We have not found blue-shifted features in the wings of 
the brightest [N 11], [O 11] and [O ill] lines, indicating that 
the faint features in the blue wings of H I lines can not come 
from emission of high-velocity ionized material. 

• They are narrower (FWHM ^ 10 km s~^) than hydro- 
gen lines (FWHM ~ 20 km s"^). FWHM has been derived 
from gaussian fits, after correcting from the underlying blue 
wing of the corresponding Balmer line, and after quadratic 
subtraction of the instrumental point-spread function. Al- 
though the relatively low velocity resolution of our spec- 
tra, which is not enough to derive precisely the value of 
the thermal width of the deuterium Balmer lines (colons in 
the FWHM indicates high uncertainties, and the low values 
of the FWHM of the deuterium lines are because deuterium 
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Figure 5. Wings of Ha to H(5 in S311. The lines are centred at 
km s~^ velocity. The dotted line of the left correspond to the 
average wavelength adopted for the D I lines. 



Table 12. Deuterium Balmer line characteristics in S 311. 

Line D I Isotopic FWHMDl FWHMHl Dl/Hl ratio 
shift (km s-i) (km s^i) (km s^i) (xlO^*) 



a -83.6 8; 20 7.0 ± 0.7 

/3 -81.4 <4: 20 8.5 ± 1.6 

7 -80.1 ~0: 20 9.4 ± 2.1 

5 -85.5 17: 21 12.9 ± 3.9 



line widths are of the order of the instrumental one) , it is suf- 
ficient to compare qualitatively with the value of the width 
of hydrogen Balmer lines (see Table lT^ . This result supports 
the idea that deuterium lines arise from a cold material with 
smaller therm al velocity, probably t he photon dominated re- 
gion or PDR (iHebrard et al.ll2nnnah. 

The detection and identification of the deuterium 
Balmer lines i n a n H ll region were first reported by 
iHcbrar d et alJ l|200 0a') from VLT/UVES data; they con- 
firmed the detection and identification of deuterium Balmer 
lines (up to Dr?) in the Orion nebula. Subsequently, 
iHebrard et ail ^l2000b^ published the detection and identifi- 
cation of at least Da and D/3 in four additional H ll regions 
(M8, M16, M20 and DEM S 103 in SMC), and confirmed 
fiuorescence as the main excitation mechanism of Dl, re- 
combination being negligible. The Di/Hi ratios presented 
in Table [T^ correspond to the intensity ratios and are upper 
limits to the abundance ratio because in addition to recom- 
bination the D I line intensities include the fiuorescence con- 
tribution that is considerably larger than the recombination 
one. 

On the other hand, lO'DeU. Ferland fc Hennevl (1200 ll) 
presented observations and a model for the emission of the 



deuterium lines in Orion, and concluded that they are pro- 
duced by fiuorescent excitation of the upper energy states 
by the far-UV radiation of the ionising star. In S 311 the 
Da/Ha and D/3/H/3 ratios are somewhat larger than in the 
case of the Orion nebula. In the l i ght o f the model out- 
lined bv lO'DeU. Ferland fc Hennevl ll200ll) . and taking into 
account that the spectral types of the ionizing stars of both 
nebulae are similar, this can be due to an additional contri- 
bution of UV radiation from other nearby cooler stars or to 
a lower UV grain extinction in S 311. On the other hand, the 
comparison of the Balmer decrements of the hydrogen and 
deuterium lines observed in our spectra follow closely the 
standard fluorescence model by lO'Dell. Ferland fc Hennevl 
( 200L . see their Figure 13) for the Orion nebula. 



9 DISCUSSION 

Few optical spectrophotometric studies of S 311 have 
been published in the literatu re. The most com- 
plete are those presented by te^le^l l|l978|) and 
IPeimbert. Torres-Peimbert fc Ravel il978ll . Unfortunately, 
these papers do not study the same slit position as ourselves 
(see §|21l, and the ionic abundances they have obtained are 
also very different. In fact, the ratio 0'''/0^''', which is an 
indication of the ionization degree of the nebula is very 
similar in these works ranging from 0.89 to 1.00; on the 
contrary, our 0^/0^^=2.8 (for t^=0.00) is much larger, 
pointing out that our slit position is nearer the ionization 
front of the nebula, where the 0'''/0'''^ ratio increases 
significantly. In fact, as it can be seen in Figure our slit 
position is on the brightest part of the nebula, coinciding 
with a filament or an ionization front. 

On the other hand, we can make comparisons with spec- 
trophotometric data of iShaver et al. ( 1983) . Their slit po- 
sition 2 is closer to ours. These authors only derived ionic 
abundances for 5 species: He^, O^"*", N+, S"*" and S^"*". They 
assumed temperatures much lower than ours, about 800 K 
lower for the high ionization species, and more than 1500 
K lower for low ionization ones. This implies, in general, an 
overestimation of the abundances with respect to us, except 
for O"*"^, which woul d be underestimate d. We have derived 
the abundances from lShaver et all (Il983l) line intensities us- 
ing the atomic data listed in Table 2] and we have obtained 
a good agreement between them and our results, except for 
S^'''/H"'" which is 0.21 dex lower t han our derived abun- 
dance. Although lShaver et alj lll983ll did not quote errors in 
the line intensities, this fact is possibly due to uncertainties 
in the measured fiux of the [S lll]A6312 line. Taking into 
account the faintness of this line in the spectrum of S 311 
(1% of H/3), we have estimated their error in the fiux mea- 
surement of [S lll]A6312 in about 50%, which is enough to 
explain the difference between the derived abundances. For 
the He"'"/H^ ratio the highest temperature assumed and to 
consider temperature fluctuations mak e us to derive a valu e 
0.13 dex higher than the one derived bv lShaver et alJ il983l) . 

Obtaining deep good-quality spectra of H ll regions lo- 
cated outside the solar circle is of great importance for de- 
riving radial abundance gradients in the Galactic disk. H ll 
regions in this part of t he Galaxy are scarce and usually 
very faint iRussei 1 120031) ■ In fact, part of the abundance 
data presented in this paper (C, N and O abundances) have 
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been included in recent works devote d to the calculation 
and modellinR of ab undance gradients jEsteban et aljlioosl : 
Carigi et al.ll2005b^ . 

In Table fTTI we compare S 311 and Orion nebula gas 
abun dances and the solar values. For the Sun: He comes 
ChriatenagikDalaEaaJdl and t he rest of elements 

Asplund. Grevesse fc SauvaJ l^200!^^ . As expected, the 



from 
from 

total abundances of S 311 are somewhat lower than those 
of the Orion nebula and the Sun because of the existence of 
Galactic radial abundance gradients and the larger Galac- 
tocentric distance of S 311. Interestingly, the Fe abundance 
is very different in both nebulae, this could be due to their 
different dust-depletion factors. 

It is important to appreciate that nebulae are 3-D ob- 
jects and the spectrum really corresponds to the integral of 
the emission contained in the column of gas covered by the 
slit area. This implies that the emission comes from a range 
of densities, degrees of ionization, temperatures and even 
extinctions within the column. However, our observations 
are limited to a small area covering a bright rim, probably 
coincident with a filament or an ionization front, precisely 
the brightest part of the nebula, where we expect the 3-D 
effects should be minimum. Nevertheless, it woul d be inter- 
esting to ma ke realistic 3-D jErcolano et al J'200 or pseudo 
3-D models ijMorisset. Stasihsk^T^Pefiaf l^OflS') for Galactic 
H II regions, combined with medium-high resolution long- 
slit spectroscopic data and narrow-band imaging in difi'erent 
emission line filters. On the other hand, 3-D effects should 
be much more severe in the case of extragalactic HII re- 
gions, where a small slit area covers a enormous volume of 
gas (several orders of magnitude larger than in G alactic H ll 
regions) . In a very recent paper, IPilvueinl ^2005^ has found 
that our data for S 311 do not fit his strong-line diagnos- 
tic diagrams for the empirical derivation of chemical abun- 
dances. This deviation is clearly due to the fact that the 
line fluxes accepted by the slit are not representative of the 
nebula as a whole. This effect has to be taken into account 
when spatially resolved observations of small zones of a neb- 
ula are used to apply empirical methods for the derivation 
of abundances. 



10 CONCLUSIONS 

We present echelle spectroscopy in the 3100-10400 A range 
of the brightest zone of the Galactic H ll region S 311. We 
have measured the intensity of 263 emission lines. This is the 
deepest spectra ever taken for a Galactic H ll region located 
outside the solar circle. 

We have derived the physical conditions of S 311 mak- 
ing use of several line intensities and continuum ratios. The 
chemical abundances have been derived using CELs for a 
large number of ions. We have determined also, for the first 
time in S 311, the C''"'' and O"*"^ abundances from RLs. We 
have obtained an average i^=0.038 ± 0.007 both by com- 
paring the O^""" ionic abundance derived from RLs to those 
derived from CELs and by applying a chi-squared method 
which minimizes the dispersion of He^/H"'' ratio from in- 
dividual lines. The adopted average value has been used to 
derive the abundances determined from CELs. 

We have compared Te(He i) vs. Te(Hi) for S 311 and 
other Galactic H ll regions finding that the paradigm 



is consistent with the observed behaviour of most of the 
objects, in contrast with what has been found for PNe. This 
result suggests that the process or processes that produce 
the abundance discrepancies in both kinds of objects could 
be different. 

We have detected four deuterium Balmer emission lines 
in S 311. These are the first detections of these lines in 
this object. Comparison with previous observations and with 
models lead us to support that fluorescence seems to be the 
most probable excitation mechanism of these lines. 
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